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We consider a many-body system of pseudo-spin- 1/2 bosons with spin-orbit interactions, which 
couple the momentum and the internal pseudo-spin degree of freedom created by spatially varying 
laser fields. The corresponding single-particle spectrum is generally anisotropic and contains two 
degenerate minima at finite momenta. At low temperatures, the many-body system condenses into 
these minima generating a new type of entangled Bose-Einstein condensate. We show that in the 
presence of weak density-density interactions the many-body ground state is characterized by a 
twofold degeneracy. The corresponding many-body wave function describes a condensate of "left-" 
and "right-moving" bosons. By fine-tuning the parameters of the laser field, one can obtain a bosonic 
version of the spin-orbit coupled Rashba model. In this symmetric case, the degeneracy of the ground 
state is very large, which may lead to phases with nontrivial topological properties. We argue that 
the predicted new type of Bose-Einstein condensates can be observed experimentally via time-of- 
flight imaging, which will show characteristic multipeak structures in momentum distribution. 
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I. INTRODUCTION 

Bose-Einstein condensation is an old and thoroughly 
studied quantum phenomenon, where a many-body sys- 
tem of bosons undergoes a phase transition in which a 
single-particle state becomes macroscopically occupied. 
This phenomenon has been observed in condensed matter 
systems and more recently in experiments on cold atomic 
gasesji*2 which provided a unique avenue to visualize the 
formation of the Bose-Einstein condensate (BEC). Bose- 
Einstein condensation is a phase transition driven mostly 
by the statistics of the underlying bosonic excitations 
and not by interactions. The statistics of basic particles 
are determined by the particle spin via the fundamental 
Pauli spin-statistics theorem;^. The spin must be integer 
for bosons and half-integer for fcrmions. 

In this paper, we discuss a cold atomic system of multi- 
level bosons moving in the presence of spatially-varying 
laser fields, which give rise to an emergent pseudo-spin- 
1/2 degree of freedom for the bosons. We emphasize 
from the outset that the symmetry operations in the 
pseudo-spin space are not related to real-space rota- 
tions and thus there is no contradiction between the 
existence of the pseudo-spin- 1/2 bosons and the funda- 
mental Pauli theorem. To "create" the pseudo-spin- 1/2 
bosons, one can use the experimental setup, proposed 
in Refs. 35.6.7.8.9 10, in which three degenerate hyper- 
fine ground states |1), |2), |3) are coupled to an excited 
state |0) by spatially varying laser fields. This "tripod 
scheme" leads to the appearance of a pair of degener- 
ate dark states, spanning a subspace which is well sepa- 
rated in energy from two nondegenerate bright states. 
The coupling between the dark and the bright states 
is very weak and will be neglected (adiabatic approxi- 
mation). The parameter which labels the dark states 
plays the role of a pseudo-spin index. This emergent 
pseudo-spin degree of freedom is similar to that stud- 
ied recently in the context of spinor condensates! n d 2 In 



particular, various aspects of the pseudo spin-1/2 boson 
physics were addresse d 13 ' 14 i 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 using two 
hyperfine states to support the internal degree of free- 
dom associated with the pseudo spin. The key distinctive 
feature of the systems studied in this paper is that the 
single-particle Hamiltonian projected onto the sub-space 
of the degenerate dark states generally possesses a non- 
Abelian gauge structure. I.e., the kinetic term of the cf- 

fective Hamiltonian has the form H^n — (pi — A) /2m, 
where A(r) is a matrix in the pseudo-spin space. In a 
recent Letter^ we pointed out that under certain con- 
ditions this non-Abelian gauge structure is equivalent to 
a spin-orbit interaction. To understand the nature of 
this interaction, we note that the dark states are eigen- 
states of an atom at rest. Once the atom moves in the 
spatially modulated laser field, the dark state label, i.e., 
the pseudo-spin index, is not a good quantum number 
and the pseudo-spin starts to precess about the direc- 
tion of the momentum. This coupling between the in- 
ternal degree of freedom associated with the dark state 
subspace and the orbital movement of the particle repre- 
sents the spin-orbit interaction. The spin-orbit coupling 
parameters can be adjusted by changing the properties of 
the spatially modulated light beams. These conclusions 
are based entirely on single-particle physics; the parti- 
cle statistics play no role. Below, we consider a many- 
particle system of bosons within this tripod scheme. Due 
to spin-orbit coupling, the degenerate ground states of 
the system correspond to non-zero momenta, leading to 
a new type of BEC, the spin-orbit coupled Bose-Einstein 
condensate (SOBEC). 

The article is organized as follows: In Section |TT] we 
introduce our model and discuss the properties of a non- 
interacting many body system of bosons with spin-orbit 
interactions. We find that, in general, the single-particle 
spectrum is characterized by two degenerate minima at 
finite momenta and we determine the transition temper- 
ature for the bosons condensing into these minima. In 



Section [TTT1 we study the effects of density-density inter- 
actions using a generalized Bogoliubov transformation 
(Subsection IIII Aj) . We show the quasiparticle excita- 
tion spectrum contains an anisotropic free particle com- 
ponent and an anisotropic sound similar to the conven- 
tional Bogoliubov phonon. By calculating the energy of 
the condensate, we find that for a system of N bosons 
the (JV+l)-fold degeneracy of the non- interacting ground 
state is reduced by the interactions to a two-fold degen- 
eracy corresponding to "left-" or "right-moving" parti- 
cles. The corresponding many-body wavefunction de- 
scribes a NOON stated suggesting that future studies 
of the SOBEC state in the context of quantum entangle- 
ment and quantum interference are highly relevant. For 
completeness, we also derive the Gross-Pitaevskii equa- 
tions for the spin-orbit coupled condensate (Subsection 
IIIIBI) . Linearizing the coupled non- linear equations in 
the vicinity of a stationary solution leads to a spectrum 
of excitations that reproduces the generalized Bogoliubov 
result. A possible experimental signature of the new type 
of SOBEC is described in Section [IV] We argue that a 
SOBEC can be observed via time-of-flight imaging, which 
will show a characteristic multi-peak structure of the den- 
sity profiles. We demonstrate that such a measurement 
generates distinct outputs for "left-" and "right-moving" 
condensates and thus can be viewed as a measurement of 
a qubit. A summary of the paper along with our conclu- 
sions are presented in Section [V] 

II. SPIN-ORBIT INTERACTING 
HAMILTONIAN AND SINGLE-PARTICLE 
SPECTRUM 

We start with the following many-body Hamiltonian 
describing spin-orbit coupled bosons, 

W = ^2 bl p I ^—1 - vp x a 2 - v'p y a 3 \ b fjp , (1) 
p;o,/3 1 ia ? 

where b ap and b ap are the creation and annihilation op- 
erators for bosons in the state with momentum p and 
pseudo-spin a = f , J,, &i are the Pauli matrices in the 
pseudo-spin space, and the parameters v and v' charac- 
terize the strength and anisotropy of the spin-orbit cou- 
pling. We reiterate that this type of spin-orbit-coupled 
Hamiltonian |T]) will appear within the recently proposed 
tripod scheme^£££&i£ in which three hyperfine ground 
states of an atom |1), |2), and |3) are coupled to an ex- 
cited state |0) via spatially modulated laser fields. The 
underlying laser-atom Hamiltonian is 

3 

H a -\ = fi |0) <0| + 1°) (Ml + h- c] , (2) 

where VLq is a constant detuning to the excited state 
and the Rabi frequencies consistent with the real- 
ization of an effective spin-orbit interaction can be 
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FIG. 1: (color online): Schematic picture of the band struc- 
ture described by Eq. @ with v/v' = 2.5 for a constant value 
of p z . The inside sheet represents the A = +1 band, while the 
outside sheet corresponds to A = — 1 and has a double- well 
structure with minima at p x = im» and p y — 0. 

taken as f2i(r) = il sin 6 cos(mv a x)e lmVbV , ^(r) = 
flsm6sm(mv a x)e lmv >>y, and fi 3 (r) = Qcos9^ with fl, 9, 
v a and Vb being constants (see, e.g., Refs. |g and [Tol for 
details). Diagonalizing the atom-laser Hamiltonian (J2|) 
via a position-dependent rotation i? MQ (r), with a € {T> I 
, bi,b 2 } and /i G {0, 1, 2, 3}, generates a pair of degenerate 
dark states 

|t) = sin^e-^l^-cos^e"^), (3) 
||) = cos (9 cos ^e-* 5 " 1 1) + cos (9 sin $^-^12) -sin0|3), 

with <fr x = mv a x and S y = mvby, and two non- 
degenerate bright states The pseudo-spin-1/2 
structure emerges when the problem is projected onto the 
subspace spanned by the pair of degenerate dark states J£ 
Applying the position-dependent rotation i?^ a (r) to the 
kinetic energy term in the Hamiltonian generates a cou- 
pling of the pseudo-spin to momentum [see Eq. I|T])] 
with v = v a cosd and v' = Vb sin 2 (0/2), in the given 
parametrization. These coupling constants can be easily 
adjusted by changing the parameters v a , Vb and 9 of the 
laser fields, which provides a knob to tune the strength 
and form of the spin-orbit interaction. 

Now, we concentrate on the generic case characterized 
by anisotropic spin-orbit interactions and assume for con- 
creteness that v > v' > 0. The trap potential and the 
inter-particle interaction are initially disregarded and dis- 
cussed in the following sections. The single-particle spec- 
trum of Hamiltonian |T]) is (see Fig. |T|) : 

Ex (p) = hi + A ^ v2 pi +v ' 2 pi ( 4 ) 

where A = ±1 labels the bands. The corresponding eigen- 
functions 4>\ p (r) — e tpr U\(x P ) are spinors with compo- 
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nents 



^Ta(Xp) 



cos 2 Xp + A 2 sin 2 x P — AAsin%p 



1/2 



\/2 [cos 2 x P + A 2 sin 2 Xp ] 



1/4 



u i\(x P ) = -iX sign [cos x P ] ^t-a(Xp), 



and 
(5) 
(6) 



where x P is the azimuthal angle in the (p x , p y )-plane and 
A = v'/v < 1. The unitary matrix U a \(x P ) diagonalizes 
the Hamiltonian JT]) (where a =f, J, corresponds to the 
pseudo-spin index and A = ±1 labels the eigenstates). 
It is obvious from Eq. ((J) that the spectrum of the sin- 
gle particle problem contains two minima at A = — 1 
and momenta p y — p z — and p x = ±mv ^ (see 
Fig. [I]). Consequently, the single particle ground-state 
is double-degenerate and the most general expression for 
the corresponding wave-function is 



*dw(r) 



1 



where wl > and wr > are the fractions of "left- 
" and "right-moving" states subjected to the constraint 
wl+wr = 1, while 0l and 0r are arbitrary phases. Note 
that by left /right-moving states we mean states with non- 
zero momentum average, (p) = ^fmve x . However, the 
corresponding average velocity vanishes (V p 7i(p)) = 0, 
so that quasiparticles characterized by these non-zero 
momentum single-particle states are not actually "mov- 
ing" , as long as the laser fields generating the spin-orbit 
coupling are maintained. Note that rotations in the man- 
ifold of the double-well ground-states are distinct from 
rotations in the pseudo-spin Hilbert space, as real-space 
and pseudo-spin coordinates are mixed up by the spin- 
orbit interaction. The two-fold degeneracy of the single- 
particle ground state is preserved if the system is placed 
in a harmonic trap. For a potential Vtrap — mio 2 v 2 /2, 
we can write the Schodinger equation in momentum rep- 
resentation: The trap potential plays the role of "the 
kinetic energy" and the real kinetic term produces a 
double-well potential in momentum space, see Fig. [T] 
The tunnelling processes connect the degenerate vacua 
in momentum space 24 . However, they do not eliminate 
the double-degeneracy of the single-particle states, which 
is protected by the Kramers-like symmetry (see Section 

Eng. 

At low temperatures, the many-body Bose system ([1]) 
condenses into the single-particle states corresponding to 
the double-well minima. The transition temperature of 
this double-well SOBEC can be calculated using stan- 
dard text-book procedures^ Let us assume that near 
and below the transition the band with A = +1 docs 
not contribute and that we can expand the spectrum 
near the minima of the band (U|) . We define the momen- 
tum q in the vicinity of the left /right minima as follows: 
p = ±mve x + q, with q <C mv. Eq. JU) leads to the 




FIG. 2: (color online): Schematic picture of the band struc- 
ture described by Eq. Q for the isotropic Rashba-type case 
with v — v' for p z = 0. The inside sheet represents the A = +1 
band, while the outside sheet corresponds to A = — 1 and has 



(7) minima a one-dimensional circle *Jv\ + Py = 



anisotropic spectrum: 



SE(q) 



2m 



The transition temperature is 



T c = 



We see that if n 1 / 3 
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(9) 



1/6 



1 — (v'/v) -C mv, our approxi- 
mation is justified and, in particular, the density of par- 
ticles in the upper band A = +1 is exponentially small. 

In the isotropic limit A = v'/v — > 1, the transition 
temperature formally vanishes. Note that in the isotropic 
case v — v' the spin-orbit term of the Hamiltonian ([T]) 
is equivalent to the Rashba models and can be reduced 
to the latter via the rotation exp (i7rc72/4) in the pseudo- 
spin space. In this case, the spectrum has minima 

mv (see Fig. [2]). 



on a one-dimensional circle 



Px 



- Pil 



The single-particle ground-state is infinitely degenerate 
and the most general expression for the corresponding 
wave- function is 

$ ring (r) = J ^y/y^ U_( x ) e l ^e^ mv ^ cosx+ysin ^\ 
o 

(10) 

where w(x) > is the angle-dependent weight of the 
Bose-condensate on a circle [J dx/(2ir)w(x) — 1] and 
<j)(x) is the angle-dependent phase. An especially in- 
teresting class of ground states corresponds to w{\) n °t 
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vanishing anywhere on the circle. In this case, the phase 
4>{x) must satisfy the constraint cf)(x + 2tt) — 4>{x) — 2irn, 
with n G Z = 7Ti(S' 1 ) being an integer winding num- 
ber. Therefore, there may exist a number of topologi- 
cally distinct ground states (characterized by the wind- 
ing number) , which can not be deformed into one another 
via any continuous transformation. We note here that a 
transition into the ring SOBEC is similar to a "weak- 
crystallization transition" discussed by Brazovsky2£ (see 
also, Refs. 28). In this case, the phase volume of fluctua- 
tions is very large, which drives the (classical) transition 
first order. Even though the transition temperature into 
the ring SOBEC vanishes in the thermodynamic limit, in 
a finite trapped system, the energy scale for the crossover 
into this state will be non-zero. 29 



III. EFFECTS OF DENSITY-DENSITY 
INTERACTION 

The most general ground-state many-body wave- 
function of a non-interacting "double well BEC" is 



N 



- y/n\(N-n)\ 



B 



iy (*) 



N- 



where n and N — n are the numbers of "left-" and 
"right-movers," B^ R arc the corresponding creation 
operators, and c„ are arbitrary coefficients satisfying 
Yin l c «| 2 = 1- I 11 the absence of spin-orbit interac- 
tion, a two-component bosonic system has a ferromag- 
netic ground-state with fully polarized pseudo-spin i 12 i 30 
We emphasize that this is not the case for the double- well 
many-body ground-state (fTTj) that describes the spin- 
orbit interacting BEC. All the arguments used for prov- 
ing the ferromagnetic nature of the ground-state for a 
two-component system^ 2 - are now irrelevant, as the real- 
space and spin components of the wave-function cannot 
be factorized due to the spin-orbit coupling. The non- 
interacting ground-state (fTTj) has an (N + l)-fold de- 
generacy. We show bellow that this large degeneracy 
is partially lifted by interactions and reduced to a two- 
fold degeneracy. We assume a density-density interaction 
H int = I J d 3 rd 3 r'n(r)V int (r - r')rt(r'), where n(r) = 
V'/'t( r )' ! /v( r ) an d Vv( r ) is the field operator, which 
is initially defined in terms of the creation/annihilation 
operators for the original hyperfine states. First, we per- 
form the position-dependent rotation i? AtQ (r) to obtain 
the effective interaction term, which has the standard 
form 
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Hint — 7TT7 



^7 ^ M nt (q)6j ip 6 ce p +q 6j 3p ,6 / 3p,_ q , (12) 
p,p',q 



where &t is the creation operator for a state with mo- 
mentum p and pseudo-spin a in the dark state subspace. 
We need to perform a second momentum-dependent 



transformation defined by ([5]) and ([6]), which introduces 
new bosonic operators labelled by the band index A — 
±1 ©: B\ p = U{ a (p)b ap and B{ p = ^ p [/ ttA (p), where 
^qa(p) = U a \(Xp) and the summation over the spin 
index a is implied. In the limit of relatively weak in- 
teractions, Vint <C mv 2 /2 (we emphasize that the spin- 
orbit coupling strength can be tuned to be arbitrarily 
strong by adjusting the properties of laser fields), the 
upper band with A = +1 is irrelevant for the low-energy 
physics. Thus, it is convenient to express the Hamil- 
tonian in terms of left /right-moving operators, defined 
as B L / R q = B_ x T (, 



q+mv ) 



Correspondingly, we have 
Ul/r a(q) = U-i a(T(q + mv)). This leads to the fol- 
lowing interaction Hamiltonian 

1 ' 

Hint = ^ ^Vint(q ( r)Si 1 pS ( T 2 p+q^3p'S (T4 p,_ q 

p,p',q {o-i} 

xC/t iQ (p)[/ Qa2 (p + q)C/] 3Q ,(p')C/ ct v 4 (p' - q). (13) 

where the prime sign in the sum over the left and right 
indices Oi — L/R = =F is restricted by the condition 
ci + o"3 = 02 + 04, i.e., the numbers of left- and right- 
movers are conserved, and q CT = q (o\ — 02)mve x . We 
stress that equation (|13p is valid in the limit of weak in- 
teractions (relative to the spin-orbit coupling) and low 
temperatures, when only single particle states with mo- 
menta in the vicinity of the two minima are occupied. 



A. Generalized Bogoliubov transformation 



V 2 
1 N 



Next, we introduce the projection operators Vn,u = 
that select the subspace characterized by n 
left-moving and (N — n) right-moving quasiparti- 
cles. The Hamiltonian can be expressed as H = 

Yn=Q^ > N.n'H'PN.n = En=o^«' An important observa- 
tion is that the Hamiltonian containing the interaction 
term (|13[) preserves the number of left- and right-movers 
and thus we can consider different "sectors," H n , inde- 
pendently. Our goal is to diagonalizc each term Ti n using 
a mean-field scheme and reduce the many-body Hamil- 
tonian to the form 



JV 



q.rr 



V N , 



(14) 

where £o(n) is the contribution of the (n, N — n) sec- 
tor to the condensate energy, while Q, a {n 1 c\) represents 
the spectrum of quasi-particle excitations. To obtain 
the mean- field result, we use a Bogoliubov- type approx- 
imation in which the operators corresponding to q = 
are replaced within each sector (n, N — n) by c-numbers, 



B L V™o e i0/2 and B R -> V^o - n e"^ 2 . Next, 
we notice that at low temperatures, the momenta of un- 
condensed bosons are q mv. Thus, we can expand the 
products of [/-vectors in (| 13[) in terms of the deviations 
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q from the minima of the energy bands 



0(x 2 )] via a generalized Bogoliubov-type transformation 



4(qi)t/ L (q 2 ) = ^(qi)^(q 2 )«f (15) 



with A = v'/v < 1 and corrections of order 0(qf 2 ) and 
2); respectively. Consequently, contributions to the 
mean-field Hamiltonian can be expanded in the small pa- 
rameter Xq — A 2 q 2 / (mv) 2 . In the zero-order approxima- 
tion, i.e., neglecting contributions of order Xq and higher, 
the mean-field Hamiltonian for the (n, N — n) sector is 



n 



N 



+ B c 



%4> 



Bo 



(l + $)e** 
Vl~F 



s( q) + 1 + g VI ~^ 2 e" 
VT^e 4 '* s(q) + 1-5 y " q 

B q 4- h.r. . (16) 



VT^2 
(l-<5)e-** 



where (5 = 2n/N - 1, B q = (#1^,-8^ is the 

annihilation operator in a spinor notation, s(q) = 
2(5 E(q)/ [no Vi n t (q)], and SE(q) is the anisotropic spec- 
trum (jHJ) near the minima. We now introduce new 
bosonic operators £>-. q = \/l — n/A_B iiq e~"^ 2 — 
^/NB R ,g e 1 ^ 2 and 5+, q = ^jNB L ,qe-^' 2 + 
\/T— n/NB^qe^^ 2 . The Hamiltonian becomes diag- 
onal for the B- -particles, which have the "free" spec- 
trum 6E(q), and has the standard Bogoliubov form 2 ^ for 

the -B+-particles. Introducing the new operators /3_. q = 



B_ q and /3_j 



Aq = -s(q) - 1 + J[s(q) + I] 1 - 1, we get 



= + J2{ n-(q)41, q /9-,q + «+(q)/3+, q /34 



with 



,(17) 



where £q ^ is the condensate energy 2 ^ in 
the zero-order approximation, f2_(q) = 
{q 2 + q 2 + q 2 [l - (v'/v) 2 ] } /(2m) is the anisotropic 
free particle quadratic spectrum and 51+ (q) = 

\J [fi_ (q) + nVhit (q)] 2 - ™ 2 V|nt (q) is an anisotropic 
sound similar to the conventional Bogoliubov phonon 
mode in a BEC. At this level of approximation the con- 
densate energy is n-independent (i.e., it is the same for 
any particular sector characterized by n left movers and 
(N — n) right movers) and, consequently, the degeneracy 
of the non-interacting ground state is preserved. In 
the first order approximation, the mean-field Hamilto- 
nian (fT6)) acquires sector-dependent corrections of order 
i q C 1. Following the above recipe, we introduce a set 

-J % 

of new operators B± q that diagonalize the B^Bq term 
in the Hamiltonian (fTB|) but not the other terms. Next, 
we diagonalize the full Hamiltonian [up to terms of order 



q = #-,q + XqDqB 



(18) 

+ XqFiqB + ,q + XqF 2 qB\__q 
(l-Al)-V 2 (B +i q + AqBl_q (19) 
+ 2;qClq-B_, q + XqC 2 qB f _ _ 



In the equations (fT5|) and p^|) we already anticipated 
that some of the terms are corrections of order Xq. 
The coefficients are determined by requiring that the (3- 
operators obey standard commutation relations [to or- 
der 0(xq)] and that the off-diagonal contributions to 
the Hamiltonian vanish. Assuming for simplicity that 
we have a point-like interaction, i.e., Vint(q) = Vint is 
momentum-independent for momenta in a range that is 
relevant for the problem, the groundstate energy in the 
(n, N-n) sector is 



£o(n) = 



8 



VintN 2 

2V 



(cos(4£) + 3) - A q (cos(4C) - 5)] } + 0{x\), (20) 



where cos 2 (£) = n/N. The relevant coefficient of the 
generalized Bogoliubov transformation (|18H19p has the 
form 



A, = -1-5 + 5^(4 + 



}y/s{4 + . 



32V«(4 + s) 
x [-4 - 5s + (4 + s) cos(4£)] + 0{x\). 



Explicitly evaluating (j2T))) with ^4 q given by Eq. (j2"Tj) 

shows that, at this level of approximation, the energy 
of the condensate becomes sector-dependent, £o( n ) ~ 



(21) 



£q + S^'ln), and is minimal for n — and n 
Thus, the density-density interaction reduces the large 
(N + l)-fold degeneracy of the ground state to a two-fold 
degeneracy. Consequently, in the limit of vanishing in- 
teractions Vint +0, the most general expression for the 
many-body wave-function is 



N. 



|*iv) 



(«)' 



1 1 vac), 
(22) 



where i«l/r represents the fraction of the left /right 
movers and 0l/r are arbitrary phases. Notice that 
Eq. (|22[) describes a fragmented or entangled BEC, un- 
less wlWr = 0. I.e., the many-body state (j2"2")l does not 
correspond to the condensation into one single-particle 
state. We reiterate that the left- and right-movers in the 
condensate have non-zero momentum, but zero velocity 
and do not actually move while the laser fields responsi- 
ble for the spin-orbit coupling are present. We also note 
that equation (|2"2"|) describes a so-called NOON stater&^L 
which is quantum correlated state with properties that 
can be exploited in applications such as quantum sensing 
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and quantum metrology. This suggests that the possi- 
bility of using spin-orbit coupled condensates as qubits 
deserves to be further investigated. 



B. Gross-Pitaevskii equations 

Let us consider the density-density interaction poten- 
tial as a contact pseudo-potential, Vi n t(r — r') = Vi n t<5(r — 
r'), where Vint — ^—a and a is the inter-atomic scat- 

/' ,11L m 

tering length. The full many body Hamiltonian can be 
written as 

# = E / ^^iwvtw (23) 

+ / d 3 r^(r)^(r)^(r)4(r), 

in terms of field operators ^(r) for the original hy- 
pcrfine states, \i G {0,1,2,3}. In Eq. (|23|) we used 

the notation = { £- + V^ rap + i?o-i } for the 

single particle Hamiltonian in the presence of a trap 
potential Vtrap j in addition to the spatially varying 
laser fields that interact with the atom, H a -i. In the 
adiabatic approximation, after projecting onto the dark 
state subspace, the first term in Eq. (j2"3"|) becomes 



E b l P { [p 2 / 2 ™ + Kr a p] 1 - vp x a 2 - v'p y a 3 } al3 bpp, 

where w ap and b ap are the creation and annihilation 
operators for bosons with pseudo-spin a = The 
interaction term is given by equation (|12p. Before 
writing down the Gross-Pitaevskii equations, let us 
summarize the three different representations used for 
describing the system of bosons interacting with the 
spatially modulated laser fields. 

i) Hyperfine states representation: This is the most 
straightforward way to describe the motion of the bosons 
and their interaction with the trap potential (Vtrap) and 
the laser fields (H a _i), as well as the density-density in- 
teraction (second term in Eq. (J22J) )■ The field oper- 
ator that creates a particle in the hyperfine state fi € 
{0,1,2,3} at point r is Wj(r), while the creation of a 
free-moving particle with momentum p is described by 



j d 3 r e zpl 'ipl(r) . By performing the position- 



ed 

dependent rotation i? MQ which diagonalizes the atom- 
laser Hamiltonian and projecting onto the dark states 
subspace we switch to the pseudo-spin representation. 

ii) Pseudo-spin representation (dark states represen- 
tation): This is the natural framework for describing 
the low-energy physics of the atomic system interact- 
ing with the laser field. The creation operator for free- 
moving particles with spin a £ {TiJ.} an d momentum 
p is 6j, p . We can define the corresponding field oper- 

ator as ip a {r) — J2 P e lpr b ap - Note that the field op- 
erators in the hyperfine and pseudo-spin representations 



are related via the position-dependent rotation, $ (r) = 

tJt 

E a ^Ma( r )V , a( r )- Diagonalizing the single-particle spin- 
orbit coupled Hamiltonian, H — [p 2 /2m + Vt rap ] 1 — 
vpx0^2 — v'pyds, generates a set of eigenstates described 
by the spinor eigenfunctions cf) an (r). The quantum num- 
ber a = ± can be viewed as labeling right (left) moving 
states. 

iii) Right/left moving states representation: This is the 
representation corresponding to the eigenstates of the 
spin-orbit coupled single particle Hamiltonian. In Section 
U we have shown that in the absence of a trap potential 
the single particle spectrum for the generic case v =/= v' 
is characterized by two minima at non-zero momenta. 
Here we show explicitly that the double-degeneracy of 
the single-particle states is a general property of the spin- 
orbit interacting Hamiltonian, protected by a Kramers- 
like symmetry. Let us use the following parametrization 
for the eigenfunctions: 



M L(r) 

*<™L( r ) 



(24) 



where a = ± and n is a set of quantum numbers. The 
components u" n (r) are the solutions of the following 
eigenproblem 



ho-v'py ivp x \ / ul n (r) 
-ivp x h + v'py J \ iau l an {r) 



E„ 



\ (25) 



where Hq = p 2 /2m + Vt rap is the Hamiltonian in the ab- 
sence of spin-orbit interaction. More explicitly, u" n (r) 
satisfy the following system of coupled differential equa- 
tions: 



V 2 , d 

2m ay 

■ d s 
-iov— — Y mv 

ox 

V 2 . d 

IT- + ^trap (r)-lV — - 

2m ay 

■ 9 r 

-iov— — Y mv 
ox 



E 



mv 



(t4(r)-ui(r)) = 0, 



E 



4(r) (26) 



(4(r)-4(r)) = 0. 



Taking the complex conjugate of (|26p with a — > — o we 
obtain an identical set of equations. Consequently we 
have 



,( r ) = [«L( r )]*: 



(27) 



and the corresponding energies are degenerate, E- an = 
E an = E n . Because ((f>-. r7n \(f> crn ) = 0, the two states 
are linearly independent. We conclude that the single- 
particle eigenstates of the spin-orbit coupled Hamilto- 
nian are (at least) double degenerate independent of the 
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symmetries (or lack of symmetry) of the trap poten- 
tial. Note that this double degeneracy is a consequence 
of a Kramers-like symmetry of the spin-orbit interact- 
ing Hamiltonian, which contains terms that are either 
quadratic in momentum, or linear in both momentum 
and spin. The creation operator for a left/right mov- 
ing particle described by the eigenstate (f> an is B\ n . The 
field operators in the pseudo-spin representation can be 
expressed in terms of B an operators as 

£ T (r) = ^[e™ ra «|„(r)B + „ + e^™ V n (r)B_ n ] , 



^(r) = Y, [ie^uU^B+r 



-imvx I 



M i„(r)B_J(28) 



where the terms with a = + and a = — correspond to 
the right and left moving modes, respectively. Finally, 
note that in the translation invariant case, T4 rap = 0, 
we introduced the eigenfunctions <fx P (r) — e lpr U\(x P ), 
with U a \(x P ) given by equations ([5]) and ©, and the 
corresponding creation operators, B\ p . We then defined 
the left /right movers for the low energy band A = — 1 
and small momenta q < mv as Bl/r q = -B(-i) =p(q+mv)- 
Alternatively, we can directly define the eigenfunctions 
</> CTq (r) in the left /right moving representation using the 
parametrization (|24[) . with no restriction for q. The cor- 
respondence between the two representations is given by: 
p = cr(q+mv) and A = — sign{q x +mv). This generalizes 
our definition of the left/right moving modes to arbitrary 
energy. Notice however, that a left (right) "moving" state 
from the high energy band A = +1 has in fact a positive 
(negative) momentum. 

To write the Gross-Pitaevskii equation in the pseudo- 
spin representation we use the standard procedure and 

calculate the commutator [i/j a (r),?{], where fi is the 
many-body Hamiltonian expressed in terms of pseudo- 
spin field operators. Using Eq. ([23]) and the relations 
between representations summarized above we obtain 



. d 



dy 



2m 



+ Vtrap(r) 



? ■ d „ 

1 +IV — G2 
OX 



(29) 



$p(r,t) + v M (lvM 2 + l^l 2 ) Mr,*)- 



Relation (|29[) , which is a system of two coupled non- linear 
differential equations, represents the time-dependent 
Gross-Pitaevskii equation for a spin-orbit coupled Bosc- 
Einstein condensate wave-function. Similar equations 
can be written in the left/right moving states representa- 
tion. For simplicity, we will address here only the trans- 
lation invariant case Vtrap = 0. The field operator for the 
left/right moving modes can be written in terms of the 
corresponding B a q operators as 



far) = £ C q (r)S, 



erq- 



(30) 



The non-interacting part of the Hamiltonian is diago- 
nal in terms of left /right moving operators, with eigen- 
values that depend on the momentum q only. At low- 
energies, these eigenvalues are given by the anisotropic 
spectrum 5E(q) = (ql + q 2 ) / (2m) + q 2 / (2m y ) with m y = 

m [l — (i//v) 2 ] . In general, the interacting Hamilto- 
nian is given by equation (|13|) . but in the low-energy limit 
we neglect all corrections of order x q = A 2 q 2 /(mv) 2 and 
higher coming from the momentum-dependent matrices 
f7 a<7 (q). In this limit we obtain 



d ~ . . 



(—id x — amv) z 
2m 



d 2 

2m,, 



d 2 \~ 
2m / 



+V int [\i} L \ 2 + \^ R \ 2 )^(r,t) 



(31) 



a.q 



where dj = d/dxj, j G {x,y, z}. The time-independent 
Gross-Pitaevskii equations can be obtained by looking for 

a stationary solution of the form i^ a (r, t) = ^ Oo .(r)e _v *, 
where fi is the chemical potential which determined by 

the condition N — J d 3 r ^IV^lI 2 + IV'rI 2 ^ > with N being 

the total number of bosons. We note that by linearizing 

-0 (T (r,t) with respect to the deviations from the station- 
ary solution we obtain an excitation spectrum consisting 
in two modes, fi±(q), identical with those found using 
the generalized Bogoliubov treatment. 



IV. EXPERIMENTAL SIGNATURE OF 
SPIN-ORBIT COUPLED BEC: MEASURING A 
SOBEC QUBIT 

A straightforward way to detect experimentally the 
new type of BEC would be to probe the momentum dis- 
tribution of the density of the particles via time-of-flight 
expansion. After removing the trap and the laser fields, 
the boson gas represents a system of free particles, each 
characterized by a certain momentum and a hyperfine 
state index. In a TOF experiment one determines the 
momentum distribution by measuring the particle den- 
sity at various times after the release of the boson cloud. 
The operator associated with a density measurement is 

P( r ) = 2~2 fJ ,i ; ji( r )i J ^( r ), where ^(r) is the creation op- 
erator for a particle in the hyperfine state \i positioned 
at point r. Determining the density profile involves a 
simultaneous measurement of p(r) for all the values of 
r G V corresponding to a ceratin region in space where 
the boson cloud is located. To insure formal simplicity, 
we consider a coarse-grained space, i.e., we treat r as a 
discrete variable. This is simply a technical trick and 
does not influence the final result. Our goal is to find the 
most likely spatial distributions of the particles at a given 
moment t after the release of the atoms. In the limit of 
large particle numbers, the actual measured density pro- 
files will involve only small fluctuations away from these 
"most likely" distributions. 
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For a system of N bosons, the result of the measure- 
ment is a set of eigenvalues {X^ t n r^}(rev) that label an 
eigenvector of the density operator 



n«{w>= n 



i 



/j, rev 



r// 



^t(r) "Hvac), (32) 



where the occupation numbers satisfy the constraint 
r n r/j = A, and the factors 1/ ^/ (n r/J )! insure the nor- 
malization to unity. Note that n rfl is an integer repre- 
senting the number of particles located in a certain "cell" 
r of the coarse-grained space. At time t after the release, 
the many-body state of N bosons that were initially in a 
BEC groundstate described by Eq. J22} is 



a 

e In 



(33) 



1 



Q2(r,i)$,(r) "11 



vac) 



where A/" is a normalization factor, a labels the left 
(a = L = — 1) and right (cr = R = +1) modes and 
||Wjv(0)) = H^iv)- The coefficients Q° are normalized 
so that \Q°(r,t)\ 2 = 1. The second summation in 

([33]) is over all possible spatial distributions of N particles 
and, in the continuous limit, it becomes a path integral. 
Equation (|33j) represent the expansion of the many-body 
wave- function in terms eigenstates (13"2"j) of the density op- 
erator. The probability ^{nr^}] of measuring a certain 
density profile n r/i is determined by the coefficient of the 
corresponding term. If we focus, for simplicity, on the 
case when there are only left (right) movers in (|2"2"|) . this 
probability is proportional to F] r |Q£(r, t)| 2nr '*/(n r|1 )!, 
with cr = L(R). The probability ^{nr^}] has a maxi- 
mum for njL = iV|Q^(r, i)| 2 corresponding, in the con- 
tinuous limit, to a stationary point of the path integral 
in equation (|3"3"|) . For large particle number, "P[{n r/1 }] be- 
comes sharply peaked at n° M and the actually measured 
density profiles will exhibit only relatively small devia- 
tions from the stationary profile. Therefore, at time t 
after release, the density of the boson cloud is 



(34) 



If both wr and u>l are non-zero, the result of a mea- 
surement will be either a "right moving" density profile 
[a = R in (|34p] with a probability wr, or a "left moving" 
profile [cr = L in (f3"4")) ] with a probability wj,, assuming 
that the two profiles are spatially well separated. We are 
not addressing here the interesting effects of the interfer- 
ence between left and right moving condensates. These 
effects are negligible if the left and right moving density 
profiles are spatially separated, but become important 
otherwise, e.g. at small times after the release. 

Next we determine explicitly the coefficients QZ(r, t) 
for the exactly solvable model of bosons with "Ising-type" 



spin-orbit coupling, v 7^ v' — 0, placed in a harmonic 
trap, Vt rap (r) = mu 2 r 2 /2. 10 In this case, the operators 
B\ from Eq. (j2"2"|) are creation operators for the single 
particle ground states 



V2 V %a 



(35) 



where </?o( r ) represents the groundstate wavefunction of 
the harmonic oscillator. The spinor (|35p is written in 
the dark state basis. Performing the position-dependent 
rotation R^ a [see Eq. ©], we can express the opera- 
tors B\ in terms of creation operators for particles in 
a certain hyperfine state located at point r, ^t(r), or 
their Fourier components corresponding to free moving 

particles, cJ k = J2 r e * kr V'/i( r )- The time evolution af- 
ter the release can be easily described in terms of time 
evolution for the cj^ operators, cj^ k (f) = exp(— ie^t) $ k , 
where = k 2 /(2m) is the free particle spectrum. Con- 
sequently, the many-body state 1 1 \1/ jv (*)) can be ob- 
tained by making in Eq. (|22|) the substitution B\ — > 



(r,t)= fcol (r')R; a (r')e 

A — ' L J a 
,k 



ik(r-r') e -iek*_ 



(36) 



Finally, introducing this expression of in equation 
(IM)) we obtain for the measured density profile the ex- 
pression 



p(r,t) 



N- 



[2tt(1 + t 2 )]^ 



sin 6 1 e !+ T 



■(« 



+ 



+ 



(1 — A cos f 
2 

(1 + Acosf 



e 1+ 



T 2 /■„ (Am^-m^)t ^ 



(37) 



where T = y/moj is the inverse characteristic length of the 
trap potential, r = u>t is time in units of u)^ 1 , 6> 6 [07r/2] 
and v a are tunable parameters characterizing the laser 
field, and v — v a cos 8 [see the paragraph containing Eq. 
(f3|)]. In equation ((37|) the density was normalized so that 
J d 3 r p(r,t) — N. The density profile for a "left mov- 
ing" density distribution (cr — —I) is shown in Fig. [3] for 
three different times after the release of the boson cloud. 
The parameters of the calculation are 8 = tt/3 and v a = 
6y/ui/m. Notice the three-peak structure of the density, 
corresponding to the three exponential terms in equation 
(|3T)) . The relative weights of the peaks are cos 4 (8/2) 
(large peak), I/2sin 2 (9 = 2 sin 2 (6>/2) cos 2 (6>/2) (middle 
peak) and sin 4 (8/2) (small peak) and their characteristic 
velocities are — cru a (f — cos 6*), av a cos9 and av a (l+cos8), 
respectively. The "left" and "right moving" distributions 
are symmetric with respect to reflection (see 

also Fig. [4|). Notice that the total momentum corre- 
sponding to a distribution described by equation (|37D 
vanishes. By analyzing the transformation ([3]) to the 




Y=8,e = 7t/4 



e = jt/3 



FIG. 3: (color online): Density of particles at three differ- 
ent moments, t\ = 0.4a; -1 , ti = 0.6a; -1 , and t$ = 0.8a; -1 , 
after both the trap and the laser fields are removed at 
t = 0. For clarity, the density distributions are shifted 
along the y-axis. This time-of-flight expansion corresponds 
to a many-body ground state [122[) and is obtained using the 
single-particle eigenfunctions for the exactly-solvable model 
of trapped bosons with Ising-type spin-orbit coupling (v 7^ 
v' = 0) 10 with v a = 6(oj/m) 1/2 . This "left moving" density 
distribution is measured with a probability Wl, while there 
is a wr probability to observe a "right moving" distribution 
which corresponds to a x — + —x reflection (see also Fig. [3J. 
Notice the characteristic three-peak structure. To resolve the 
BEC peaks, the spin-orbit coupling energy scale should be 
larger than the trap level spacing, i.e., mv 2 ^> lu. In the op- 
posite limit the phenomenon of real-space BEC separation is 
smeared out by finite-size effects (|22p . 



dark state basis, we observe that sin 6* is the coefficient of 
the hyperfine state |3). Consequently, the middle peak in 
the density distribution ([57)1 consists of particles in this 
particular hyperfine state. The other two peaks contain 
mixtures of states |1) and |2). A state-selective measure- 
ment of particles in the hyperfine state |3) would reveal 
a single peak structure moving to the left or to the right 
with a velocity v — v a cos 9. The dependence of the den- 
sity profile p(r, t)/N on 9 and on the ratio 7 = v a / y^oj/m 



for r = (x, 0, 0) and t 



is shown in Fig. |4j 



V. SUMMARY AND CONCLUSIONS 

To summarize, in this article we have introduced and 
discussed in detail a new type of many-body system con- 
sisting of pseudo spin-1/2 bosons with spin-orbit inter- 
actions. We have shown that at low temperatures the 
system condenses into a new type of entagled BEC, the 
spin-orbit coupled Bose-Einstein condensate (SOBEC). 
The novelty of this state stems from the coupling of an 
internal degree of freedom, the pseudo-spin created as a 
result of an atom interacting with a spatially modulated 




-18-12 -6 6 12 

x 



4 4 8 12 16 
x 



FIG. 4: (color online): Density profiles p(r,t)/N for r = 
(a;, 0,0) and t — oj~ 1 . The position x is measured in units of 
r - . Left panel: "right moving" versus "left moving" distri- 
butions. Notice that the "center of mass" of the distributions 
is always at x — 0. Middle panel: Dependence on the angle 
9. At small angles all the weight concentrates in the large 
peak which is centered near x — 0. In the limit 9 — > n/2 the 
strength of the SO interaction vanishes v — > and the present 
analysis is not valid. Left panel: Dependence on the relative 
strength of the spin-orbit coupling, 7 = v a j yw/m, To re- 
solve the peak structure, the spin-orbit coupling energy scale 
should be larger than the trap level spacing. In the opposite 
limit interference effects become important (see main text). 



laser field, to the real space motion of the particles. As 
a result, the single-particle spectrum is characterized by 
degenerate minima at finite momenta and, consequently, 
the bosons condense at low temperatures into an entan- 
gled quantum state with non-zero momentum. For an ar- 
bitrary spin-orbit coupling, the single particle spectrum 
has a double-well structure in momentum space (see Fig. 
[U with minima at non-zero momenta. In this case, a 
system of N non-interacting bosons is characterized by a 
large (N + l)-fold degeneracy of the many-body ground 
state. Weak density-density interactions reduce this large 
degeneracy to a two-fold degeneracy. The corresponding 
ground state wave-function describes a superposition of 
left-moving and right-moving condensates with weights 
wl and wr = 1 — u>£, respectively. Performing a time- 
of-flight expansion of the condensed bosons results in a 
characteristic three-peak structure (see Fig. |4]). The to- 
tal momentum of the density profile is identically zero, 
but the peaks are moving along the x-direction with ve- 
locities proportional to the k- vector of the laser field mod- 
ulation in that direction. The probability of measuring 
a left- (right-) moving condensate is Wl (wr) and the 
signature of a left- (right-) moving state consists in the 
middle and small peaks moving left (right), while the 
large peak moves in the opposite direction. 

In conclusion, the spin-orbit coupled BEC can be 
viewed as a state occurring at the interface between spin- 
tronics and cold atom physics, with nontrivial properties 
that have a significant potential for applications. We note 
here that the ground-state of the double- well SOBEC [see 
Eq. (|22]l] represents a NOON state^M! similar to those 
recently proposed for the construction of a gravimeter 
bases on atom interferometry^ Therefore, the study of 
a SOBEC state in the context of quantum entanglement 
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and quantum interference is highly relevant. In addition, 
the double degeneracy associated with the pseudo-spin 
degree of freedom makes this state a natural candidate 
for a qubit. A possible way to measure such a qubit was 
described in the last section. Time-dependent laser fields 
[similar to those, which lead to the spin-orbit-couplcd 
Hamiltonian (JTJ)] could be used as "gates" to perform uni- 
tary rotations in the space of degenerate ground states. 
Note that the coupling of the spin to the orbital motion 
yields a protecting mechanism against decoherence, due 
to momentum conservation, and suggest that the spin- 
orbit coupled condensates are interesting candidates for 



fault tolerant quantum computation. Finally, we note 
that for a symmetric Rashba-type spin-orbit coupling the 
system is characterized by a single-particle spectrum that 
has a continuous set of minima along a circle in momen- 
tum space. This results in a huge degeneracy that may 
lead to possible phases with non-trivial topological prop- 
erties, making the study of the symmetric SOBEC a po- 
tentially very interesting problem. 

V.G. is grateful to C. Clark, I. Satija, I. Spielman, and 
J. Vaishnav for useful discussions and JQI for financial 
support. 
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